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Abstract 

We introduce a scheme to include many-body screening processes explicitly into a set of self- 
consistent equations for electronic structure calculations using the Gutzwiller approximation. The 
method is illustrated by the application to a tight-binding model describing the strongly correlated 
7-Ce system. With the inclusion of the Sd-electrons into the local Gutzwiller projection subspace, 
the correct input Coulomb repulsion Uff between the 4/-electrons for 7-Ce in the calculations can 
be pushed far beyond the usual screened value U'^^J and close to the bare atomic value U^f^ ■ This 
indicates that the d-f many-body screening is the dominant contribution to the screening of Uff 
in this system. The method provides a promising way towards the ah initio Gutzwiller density 
functional theory. 
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I. INTRODUCTION 

Over the past several decades, density functional theory (DFT) with the local density 
approximation (LDA)j2| has been very successful in electronic structure and total energy 
calculations in many systems. Meanwhile, several hybrid approaches based on the combina- 
tion of LDA with many-body techniques have been proposed to overcome the limitation o " 



LDA in strongly correlated electron systems. Among these hybrid approaches, LDA+U 



m 



is the most widely used method. Using a more detailed treatment of electronic correlation 
effects, LDA plus dynamical mean field theory (LDA+DMFT) js, 6| takes the local quantum 
fluctuations into account and can calculate both ground state and excited state properties. 
However, due to the large computational load of DMFT, LDA+Gutzwiller method has re- 
cently been developed to calculate the ground state properties of correlated systems[7|,|8|. All 
practical calculations with the above hybrid approaches use as input the screened Coulomb 
repulsion U and Hund's coupling J parameters for the correlated orbitals, which have to 
be estimated in advance using experiment data, constrained-LDA calculations joj] or random 
phase approximation (RPA) 
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Different from the above hybrid approaches which require prior determination of the 
screened Coulomb repulsion U and Hund's coupling J parameters, a Gutzwiller density 
functional theory (GDFT)[l3| has recently been proposed as an ab initio approach which 
directly takes the Coulomb integrals of the local orbitals and incorporates the screening 
process explicitly through a self-consistent solution of the many-electron wave function. In 
the GDFT, a Gutzwiller form is adopted for the variational wave function to go beyond single 
Slater determinant-based approaches, while the effective single particle picture is retained 



through the Gutzwiller approximation |l4|-|l6 1 . 

In this paper, we address an important yet unresolved question in the GDFT, i.e., how 
to choose the local subspace for the Gutzwiller operator such that the most relevant many- 
body screening processes can be captured in a self-consistent calculation. An orthogonal 
tight-binding model for the strongly correlated 7-phase Ce is chosen as a prototype to 
examine the effects of including main onsite screening 5(i-channels in addition to the localized 
4/-orbitals into the Gutzwiller local subspace. In many LDA+J7, LDA+Gutzwiller and 
LDA+DMFT calculations such screening processes are usually bypassed by introducing the 
screened interaction parameters (e.g., U and J parameters) as input between the minimal 
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correlated orbitals (e.g., 4f-orbitals for Ce). Such screened U parameter is usually much 
smaller than the bare atomic value since the screening effects are assumed to be absorbed 
in it. By the explicit inclusion of some many-body screening processes into self-consistent 
calculations, one may expect to use the partially unscreened larger input U parameter for 
the calculations in order to obtain the same physical effects. And the fully unscreened 
U parameter (i.e., bare atomic value) should be used as input for the calculations if all 
the important many-body screening processes are explicitly dealt with. In this paper, we 
show that a much larger input U parameter close to bare atomic value is needed for the 
calculations with both 5d and 4/-orbitals in the Gutzwiller local subspace in order to retain 
the same physical effects as one uses the screened U parameter with minimal 4/-orbitals 
in the Gutzwiller local subspace for 7-Cc. Therefore the treatment of 5d-4f interactions 
on the Gutzwiller level is indeed a many-body screening process. Including the many-body 
screening effects in a self-consistent way is critical for a predictive first principles theory. 
Such calculations can also provide useful insights on how to determine accurate screened 
interaction parameters. 



II. METHOD 



The electron Hamiltonian for 7-phase Ce is written as 

n^no + 'Hi + 'H2 (1) 

where the bare paramagnetic band Hamiltonian is 

^^0 = y ] ^iajli^iac^iPcr + ^iaC^iaa^ioicr (2) 

tiaj/s is the electron hopping element between orbital a at site i and orbital /3 at site j. Sia 
is the orbital level, a and j3 run over all the basis set orbitals, i.e., 6s, bd and 4/-orbitals 
of Ce. c^(c) is the electron creation (annihilation) operator, a is spin index. The typical 
simplified onsite term for 4/ electrons is 

i,77'e{4/},(7(^)7^(7'^') 

where only density-density type interactions are included. Uff is the Coulomb repulsion 
energy between the localized 4/-electrons. 7 is the orbital index in the Gutzwiller local 

3 



subspace. {4/} denotes the set of 4/-orbitals. In our present model, an additional onsite 
many-body interaction term 

^2 = Ufd ^ cJ^^Q^^c|y^,Ciy^' (4) 

«,7G{4/},7'G{5d},cro-' 

is also included where U fd is the Coulomb repulsion energy between the localized 4/-electron 



and Sd-electron as introduced in the Falicov-Kimball model 
space in this model includes 4/ and Sd-orbitals. 

We use a variational wave function of the Gutzwiller form 



17|. The Gutzwiller local sub- 



G2 



with the Gutzwiller approximation to calculate the expectation values of the Hamiltonian 
\l/o is the uncorrelated wave function and G is the Gutzwiller projection operator. 



The Gutzwiller approximation renormalizes the expectation values of one-particle operators. 
The renormalization factors have analytic expressions, however, they are usually quite com- 
plicated. These expressions can be simplified by introducing a set of rotated orbitals 
in the Gutzwiller local subspace following Ref . [16|] , such that the local single particle density 
matrix is diagonalized, i.e., 

h\ta = ^U*^^'c\~t'a (6) 
7' 



\l'o ) for a general operator O. 



Here we define {Oj = O 

We can rewrite the Hamiltonian (Eq{T]) of the system in terms of the rotated local natural 
basis set orbitals 

Here we define h}^^ = c\^^ for the orbitals outside the Gutzwiller local subspace and the set 
of £ia, and Ul^y in the new representation can be obtained from their values Uajp, Eia, 
and f/^y in the original Hamiltonian through the basis transformation. 
We introduce a Gutzwiller operator in the following form 

G = e-^^^'^^^\^'^^^'\ (9) 



where \ J^i) is the Fock state generated by a set of 



7cr 



with rq^^ = {J-'i \ni^fj\ J-'i), which identifies whether there is an electron with spin a occupied 
in orbital 7 for Fock state J-" at the i^^ site. Qijr = 1 for empty and singly occupied config- 
urations because in these cases there are no electron-electron repulsion involved. According 



16|, the expectation value of the electron Hamiltonian for 7-Ce can be expressed 



to Ref. 

as 

(^)g = i^c^aZ^aP^p + eJap) (/lL,7^k/3<x)^ " J] ^7^?<x + ^^^^ ^^^^ 

where is the Fourier transformation coefficient of Uajp at crystal momentum k. The site 
indices are dropped since there is only one atom in the primitive unit cell of 7-Ce. We define 
the Gutzwiller orbital renormalization factor Zaa = 1 for orbitals outside the Gutzwiller 
local subspace. Following ref [3], 

1 

n \ 



Y.VM^\{J'\hUr)\' (12) 



for the orbitals inside the Gutzwiller local subspace. The energy of configuration {J-") is 
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70- (7(t)^(7'o-') 



^ (13) 



Here pjr is the occupation probability of configuration IJ-"), which satisfies the following 
constraints 

In practical calculations Hq is obtained by downfolding the first-principles LDA band struc- 
ture. Therefore some contributions to the total energy of the system from "Hi and have 
already been taken into account by "Ho in a mean-field way. Such contributions are com- 
monly referred to as the double counting term which should be subtracted in the expression 
for the total energy. Therefore the total energy per unit cell of the system is given by 

Et = {n)a - Ed.c. (16) 



Following the treatment in LDA+U calculations 18|, ll9| we choose Ed c. to be: 



Ed.c. = ^UffNf {Nf - 1) + Uf.NfN, 



(17) 



where Nf{Nd) is the total number of 4f{5d) electrons. 

Minimization of the total energy Et with respect to the band wave function tl^nko- and the 
local configuration occupation probability pjr under the set of constraints given by EqjMl 
and Eq|T5] yields the following set of equations which need to be solved self-consistently, 



'^eff IV'nkcr) — Cnkcr I'^nkcr) 



where 



(18) 
(19) 

(20) 



7,(T 



and 



M 



7. 2JnO^(l-ny V 



(21) 



The effective single-particle Hamiltonian Eqj20] has been shown to describe the Landau- 



Gutzwiller quasiparticle bands 



20|, and the square of the 2;-factor corresponds to the quasi- 



particle weight [Sl]. The local orbital chemical potentials are given by 



dz. 



/^7o 



7(T 



Uf.Nfh 



/d^V^[7e{5d}] 



with 



e7<^ = Y {^P^^P (^hl^^hk^^) + c.c. 

k,/3 



(22) 



(23) 



rjaa is the Lagrange multiplier associated with the constraint of EqJTH I[x] is an indicator 
function which equals 1 if x is true and otherwise. 

To get the self-consistent solution of Eq{T8] and [19] with the constraints of EqJTl] and [T5| 
one starts with some initial guess of {z^a} and {fi-ya}- The band wave functions {ifjuka} can be 
obtained straightforwardly by diagonalizing the effective single electron Hamiltonian "H^ff- 



r xw L r xr xw l r x 



FIG. 1: (Color online) Upper left panel shows the comparison between ab initio plane-wave based 
band structure (red thick lines) and the downfolded tight-binding result (black thin lines). Other 
three panels show projected weight of 6s, 5d and 4/-orbitals for each band. 



Then Eq{T9] with constraints of Eq{Tl] and [15] can be solved by adjusting {r/^o-} such that 
the real symmetric matrix Ai generates the lowest-lying eigenvector which satisfies EqJT5l 
{z^o-} and {^■ya} would be updated according to the solution of EqUH] until convergence is 



reached 



2l|. 



III. RESULTS AND DISCUSSIONS 

The bare band Hamiltonian "Hq of 7-Ce is obtained by downfolding the DFT-LDA band 
structure from the converged large basis set to a minimal basis-set (6s, 5d and 4/) tight- 
binding representation with the recently developed QUasi- Atomic Minimal Basis-set Orbitals 



(QUAMBOs) scheme 22l425| . FigJUshows that up to 2eV above Fermi level the tight binding 



band structure agrees very well with the LDA result from VASP[26|] calculations. The orbital 
projected weight for each band has also been plotted, confirming that 4/-contributions are 
dominant near the Fermi level. The total number of 4/ electrons is found to be 0.93, which 
is fairly close to one. 

In order to single out the many-body screening effects of the 5(i-channels, we first study 
the system with 7^2 treated at the Hartree-Fock mean-field level. The double counting term 
in EqJT7]represents the Hartree-Fock mean field interactions between the 4/ and Sd-electrons 
contained in Tio, therefore, it exactly cancels the contribution from the ^2 at the Hartree- 
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Fock level and the whole electron Hamiltonian of the system reduces to the usual form with 
many-body Coulomb interactions only between the 4/-electrons (i.e., Hi). Proceeding with 
the Gutzwiller solution of "Hi, we include in our calculations 4/-shell occupations of /°, 
and resulting in 106 local /-configurations. Fig |2] shows the variation of orbital renormal- 
ization factors and orbital occupations with increasing Coulomb repulsion parameter Uff. 
Due to cubic symmetry of 7-Ce, the 4/-orbitals split into three groups with degeneracy of 3, 
3, and 1. It can be seen that all the 2;-factors initially decrease with increasing Uff. A transi- 
tion occurs near Ujj = 6eV, beyond which the z-factors of the two sets of 3-fold degenerate 
4/-orbitals increase sharply while the z-factor of the singly degenerate 4/-orbital decreases 
rapidly. Meanwhile, the orbital occupations exhibit similar variations near Ujjr = 6eV. 
The occupations of the triply degenerate orbital groups vanish swiftly as f/// exceeds 6eV, 
while the singly degenerate orbital quickly approaches half-filling. These results support 
the physical picture of orbital-selective Mott transition 27|, |28|] in which the /-electron in 
the system is redistributed among the 4/-orbitals such that the singly degenerate 4/-orbital 
approaches Mott localization with a half-filled band, while the remaining 4/-orbitals become 
empty due to strong /-electron correlation effect. However, it has also been shown that a 
finite hybridization between the narrow correlated /-bands and the wide conduction bands 
(mainly 5d bands in this system) can suppress the Mott transition by Kondo screening 29| . 
The screened UfJ for 7-Ce has been calculated using constrained-LDA method and yields a 
value of 6eV which interestingly coincides with Uj^ 30|] . Therefore we choose the transition 
point in the z-U curve as a fingerprint to identify the correct input Uff parameter for 7-Ce 
in the following calculations of different levels. We should emphasize that LDA itself can 
not capture the many-body Kondo screening effect. The constrained-LDA method gives 
reasonable estimation of the screened UjJ due to other processes, e.g., dielectric screening. 

We now study the effect of the many-body screening when the 5(i-electron screening 
is explicitly included by treating the ^2 and Tii on equal footings with Gutzwiller varia- 
tional wave functions constructed on the Gutzwiller local subspace consisting of both 4/ 
and 5(i-orbitals. Any deviation of the calculation results from the mean field calculation de- 
scribed above will reflect the effect of many-body correlation effects arising from including 
5(i-electrons in the Gutzwiller projector. Since we do not expect any significant correlation 
effect for 5(i-orbitals, we include all the local 5(i-configurations into account, which pushes 
the dimension of the local subspace to 108544. By taking advantage of sparse matrix tech- 
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FIG. 2: (Color online) Orbital renormalization factors Zr^ (left axis, blue solid symbols) and orbital 
occupations (right axis, red open symbols) as functions of Coulomb repulsion parameter between 
4/-electrons Uff. Square and circle are for the first two groups of 4/-orbitals with degeneracy of 
3. The third singly degenerate group is symbolized by up triangle. For paramagnetic solution we 
have = z^^ = and = n[^^ = n^^. 

niques, a typical self-consistent calculation can still be done quite fast on a single processor. 
A set of Coulomb parameter Ufd values of {leV, 3eV, 5eV} are chosen with the successive 
inclusion of different number of many-body screening 5(i-channels to investigate the details 
of the many-body screening process. Our calculation results are shown in Figj3] for the or- 
bital renormalization factors z-^ and in Fig JH for the orbital occupations One can see that 
the curves shift towards larger values of Uff in a nonlinear fashion with increasing Coulomb 
parameter U fd and increasing number of many-body 5(i-screening channels, indicating that 
the many-body screening is not a linear process. The correct input Uff parameter for 7-Ce 
in the calculations which corresponds to the transition point of the z-U curve is pushed 
much closer to the bare Coulomb repulsion energy between 4/-electrons (~23eV based on 
the neutral atomic 4/-orbitals) with Ufd ~ ^eV and all the Sd-channels participating in 
many-body screening process. It should be pointed out that we encounter some numerical 
problem in searching solutions of the system on the smaller Uff side, however, this numerical 
difficulty will not affect our conclusion. 

To shed some light on the many-body screening mechanism, it is instructive to examine 
the variations of the local electronic configuration occupation probabilities {pjr}, which are 
the new variational degrees of freedom in the Gutzwiller treatment. However, a straightfor- 
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FIG. 3: (Color online) The dependence of the orbital renormalization factor on the Coulomb 
repulsion parameter between the 4/-electrons Uff, that between the 4/ and 5d-electrons Ufd, and 
the number of many-body screening 5d-channels explicitly included in the Hamiltonian, which is 
1 for the left panel (one of the eg-orbitals), 2 for the middle (two e^-orbitals) , and 5 for the right 
(two e^-orbitals and three t2g-orbitals). The set of values of Ufd scanned in the calculation are leV 
indicated by black solid symbols with solid lines, 3eV by red half-filled symbols with short-dashes, 
and 5eV by blue open symbols with short dots. The three distinct z-factors of 4/-orbitals are 
depicted by the same symbols as in Figl2]with the average 2;-factor for the participating d-channels 
by down triangles. 
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FIG. 4: (Color online) The dependence of the orbital occupation on the Coulomb repulsion 
parameter between the 4/-electrons Uff, that between the 4/ and 5d-electrons U fd, and the number 
of many-body screening Sd-channels explicitly included in the Hamiltonian. Plot settings are same 
as FigiSl 



ward analysis of the local electronic configurations can be rather tedious since its dimension is 
usually very large. In Figj5]we present the sum of the occupation probabilities of the grouped 
local configurations according to the number of /-electrons occupied, which is closely related 
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FIG. 5: (Color online) The sum of occupation probabilities of some group of local electronic 
configurations as a function oiUjf and Ufd with all the 5d and 4/-orbitals in the Gutzwiller local 
subspace. Squares are for the configurations with no /-electron occupied, circles for those with 
one /-electron, and up triangles for those with two /-electrons. [/j^=OeV is indicated by black 
solid symbols with solid lines, 3eV by red half-filled symbols with short-dashes, and 5eV by blue 
open symbols with short dots. The results with Ufci=0 are same as those with only 4/ orbitals 
included in the Gutzwiller local subspace since there are no Coulomb interactions between 5d and 
4/-electrons. 

to the z-factors for the 4/-orbitals. The local configurations with two /-electrons occupied 
can survive to larger Uff and the local configurations with one /-electron occupied will ap- 
proach a full occupation at higher Uff hj increasing Ufd- This tendency is consistent with 
the behavior of the 2;-factors for the 4/-electrons. Another way to look at the many-body 
screening mechanism is to investigate the response of the screening 5(i-electrons when the 
many-body screening is turned on. Figj3] shows that the z-factors for Srf-orbitals are smaller 
than one, but are still quite large (> 0.8). It implies that the Sd-electrons rearrange them- 
selves in the local subspace in a way deviating from Hartree-Fock mean field solution in 
order to effectively screen the interactions between the 4/-electrons. 

By including Ce 5d and 4/-orbitals into the Gutzwiller local subspace for the construction 
of Gutzwiller wave function., i.e., having 5(i-electrons treated in the same Gutzwiller level 
as 4/-electrons beyond static Hartree-Fock mean-field approximation, we observe that the 
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equivalent Ufj is pushed towards its bare atomic value in 7-Ce. Therefore local many body 
response in the 5rf-channels is the dominant contribution in the screening process for 4/- 
electrons in bulk Ce. This justifies our method in using single-site Gutzwiller treatment 
in this system. In some cases where significant screening effect comes from the atomic 
environment [31j . we expect th e g eneralized cluster-Gutzwiller treatment in parallel with 



cluster-DMFT would be useful 



32|. 



IV. CONCLUSION 



A method for explicitly including many-body screening process into self-consistent cal- 
culations based on Gutzwiller variational wave function and Gutzwiller approximation has 
been proposed and applied to the tight-binding electronic structure calculation of 7-Ce. The 
resultant Coulomb repulsion parameter between 4/-electrons Uff can be much larger than 
the corresponding screened Uff with moderate value of Coulomb repulsion parameter be- 
tween 4/ and 5(i-electrons Ufd- Assuming the dominant many-body screening effect for the 
strongly correlated orbitals (e.g., 4/-orbitals for Ce) can be explicitly taken into account in 
this local approach, it is therefore of great interest to examine how the current scheme per- 
forms in full DFT-based self-consistent calculations where all the Coulomb repulsion energies 
between various orbitals in the Gutzwiller local subspace are directly obtained from atomic 
orbital integrations. Since the atomic orbitals are generally non-orthogonal, it is also highly 
worthwhile to extend the scheme from orthogonal orbitals to non-orthogonal orbitals. The 
proposed scheme here furnishes a promising way towards the realization of a self-consistent 
ab initio GDFT. 
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